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The methods used till now to calculate the neutrinoless double beta decay matrix elements are: the 
Quasiparticle Random Phase Approximation (QRPA), the Shell Model (SM), the angular momen- 
tum projected Hartee-Fock-Bogoliubov approach (HFB) and the Interacting Boson Model (IBM). 
The different approaches are compared specifically concerning the the angular momenta and parities 
of the neutron pairs, which are changed into two protons by the 0^/?/3-decay. The QRPA and SM 
involve about the same angular momentum and parity neutron pairs, while the HFB is restricted 
to 0+, 2+, 4+, . . . , and the IBM to 0+ and 2+ nucleon pairs. 

PACS numbers: 21.60.-n, 21.60.Jz, 23.40.-s, 23.40.Hc, 



I. INTRODUCTION 



Neutrinoless double beta decay (Ov/3/3 decay) allows to determine whether the neutrino is a Majorana or Dirac 
particle and gives also a value for the absolute scale of the neutrino masses and information on coupling constants 
and masses of physics beyond the Standard Model for Grand Unified Theories (GUT's), Supersymmetry (SUSY) 
and models with extra dimensions (see e.g. [2-13|)- The different CV/3/3-decay transition matrix elements are as 
important as the data on fV/3/3 decay to determine the neutrino mass and other parameters beyond the Standard 
Model. In the last three years the reliability of the transition matrix elements has been greatly improved (see e.g. 
Pt3 )■ The groups in Tubingen, Bratislava and Jyvaskyla are using mainly the Quasiparticle Random Phase 
approximation (QRPA), while Poves and collaborators [8l IT^Tlll| use the Shell Model (SM), Tomoda, Faessler, Schmid 
' and Gruemmer [l4j and Rath and coworkers jT^-[l3 use the angular momentum projected Hartree-Fock-Bogoliubov 
method (HFB), and Barea and Iachello [l8| use the Interacting Boson Model (IBM). 

The QRPA [{J 0] has the advantage to allow to use a large single-particle basis. Thus, one is able to include to 
each orbit in the QRPA model space also the spin-orbit partner, which guarantees that the Ikeda sum rule [19( is 
fulfilled, that is essential to describe correctly the Gamow- Teller strength. The SM Q is presently still restricted to a 
nuclear basis of four to five single-particle levels for the description of double beta decay. Therefore, not all spin-orbit 
partners can be included and, as a result, the Ikeda sum rule is violated. On the other side, the shell model can in 
principle take into account all many-body configurations in a given single-particle basis. The excited states in the 
QRPA for spherical even-mass nuclei include 'seniority' (the number of broken quasiparticles) states two, six, ten, 
. . . and for the ground state correlations 'seniority' zero, four, eight, . . . configurations. The SM takes for the ground 
state seniority zero, four, six, eight, and for the excited states seniority two, four, six, ...into account. Since 
the SM uses the closure relation for calculating the 0i//3j3 matrix elements, one needs only the configurations for the 
initial and final even-even nuclei in the + ground states. In contrast to the QRPA, the SM includes also states of 
the seniority six, but on the other side the SM violates the Ikeda sum rule due to the small single-particle basis. The 
effects of the seniority six configurations and of the violation of the Ikeda sum rule are also studied in this work. 

The difference of the HFB and the IBM approach to the QRPA and the SM treatment is for the later two the 
higher flexibility to include more angular momenta and parities for the two neutron pairs, which are changed into two 
protons in the Qv(3(3 transitions. The HFB approach with angular momentum projection after variation allows only 
neutron pairs with J 71 " = + , 2 + , 4 + , . . . or in general (— ) J = 1, positive parity states for real coefficient and no parity 
mixing in the Quasiparticle transformation [20L I211J (see equation ([TS])). The IBM is restricted to neutron "S" and 
"D" pairs, i. e. two neutron states with angular momenta + and 2 + , which is a restriction not found in the more 
flexible QRPA and SM approaches as shown in this work. 
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II. THE DIFFERENT MANY BODY APPROACHES AND THE TRANSFORMED NEUTRON PAIRS 

Here a short characterization to the different many body approaches: QRPA 0-0|, SM H,[l(j, HFB [IE EE E3, IH 
and the IBM p]| is given. 

In the QRPA one starts from the transformation to Bogoliubov quasiparticles : 



a] = met - Via. 



(1) 



The creation cj and annihilation operators of time reversed single-particle states cj are usually defined with respect 
to oscillator wave functions [EJ. The single-particle energies are calculated with a Woods Saxon potential [f|. The 
single-particle basis can include up to 23 nucleon levels (all single-particle states without a core up to the ii 3 / 2 level) 
for the protons and also for the neutrons. But the QRPA results for the 0vf3/3 matrix elements turn out to be stable 
as a function of the basis size already for smaller basis sets (from 6 or 7 levels and larger, respectively). 

The excited states \m) with angular momentum J in the intermediate odd-odd mass nucleus are created from the 
correlated initial and final + ground states by a proton-neutron phonon creation operator: 



\m) = Ql\0 + ); HQl\0+) = E m Ql\0+). 



(2) 



(3) 



which is defined as a linear superposition of creation operators of proton-neutron quasiparticle pairs: 

Al = [atal]j M , (4) 



[A a , Afa] = 5 a .J3 + X, 



(5) 



The "scattering" terms x in ([5]) are put to zero in the QRPA and are included approximately in Renormalized 
QRPA (RQRPA), which for the first time was employed to neutrinoless double beta decay in the PhD thesis of 
Schwieger in Tubingen 22]. The RQRPA includes the Pauli principle for the Fermion pairs and reduces the ground 
state correlations. The many body Hamilton equation ([2|) yields the algebraic QRPA equation for the the coefficients 
X™ and Y™, and the excitations energies E m . The approach with x = is the so called Quasi Boson Approximation 
(QBA). 

The inverse 0^/3/3 lifetime for the light Majorana neutrino exchange mechanism is given as the product of three 
factors, 



G°" I AT 



0i/| 2 2 



(6) 



where G°" is a calculable phase space factor, M 0v is the 0^/3/3 nuclear matrix element, and mpp is the (nucleus- 
independent) "effective Majorana neutrino mass" which, in standard notation [23| , reads 



TO/3/3 



(7) 



with mi and U e i being the neutrino masses and the v e mixing matrix elements, respectively. 

The expressions for the matrix elements M 0l/ and the corresponding OV/3/3 transition operators are given, e.g., in 
Rcf. [|: 



M (0u) = M 0u _ (£K)2 M 0, 
9A 



(8) 



M K » = (0+\H K (r 12 ,E^)\0 + ), H K (r 127 E" J l) = ^ r f K (qr 12 ) 

Kg A Jo 



hx{q 2 )qdq 



q + E^-iEi-Ef)/! 



(9) 
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a) 



b) 




FIG. 1: The ring diagrams of type a) are with the "vacuum" included in the ground state correlations in the QRPA. This 
means that the seniority, or in deformed nuclei the quasiparticle number, zero, four, eight, twelve, . . .excitations are included 
in both the QRPA and the SM. But, e.g., diagram b) is not included in the QRPA but is considered in the SM, where seniority 
4, 6, 8, 10, 12, ... are included in the ground state correlations §, ITol. HT1] . However, contributions of the seniority six and ten 
admixtures to the neutrinoless double beta decay probability are suppressed according to the philosophy of the Random Phase 
Approach, that the ring diagrams are the leading ones (see also figure 4). 

Here, "K" labels GT (Gamow Teller), F (Fermi) and T (Tensor) contribution. The functions fx (qr^) are spherical 
Bessel functions jo (qr 12) for F and GT and ^'2(9^12) for T. The function hxiq 2 ) is defined in appendix A of Simkovic 
et al. @, and are the energies of states in the intermediate nucleus. The nuclear radius is R — r^A 1 ^, and 
the parameter tq is chosen here to be tq = 1.2 fm. Thus, the matrix elements M 0v are increased by about 10 % as 
compared with our earlier calculations, where we chose = 1.1 fm. 

The SM approach has been applied by the Strasbourg-Madrid group [l(| to neutrinoless double beta decay Q 
using the closure relation with an averaged energy denominator. In this way one does not need to calculate the states 
in the odd-odd intermediate nuclei. The quality of the results depends then on the description of the + ground 
states in the initial and final nuclei of the double beta decay system, e.g. 76 Ge — » 76 Se. The Ov/3/3 transition matrix 
element © simplifies as shown in Ref. [ll[ equations (5) to (11). Since the number of many body configurations is 
increasing drastically with the single-particle basis, one is forced to restrict for mass numbers A = 76 and A = 82 in 
the SM to the single-particle basis lp3/2,0/ 5 / 2 , lpi/2 an d 0g 9 / 2 . In ref. @ the 82 Se nucleus is calculated in the SM 
for five basis single-particle levels including also 0/ 7 / 2 - For the mass region around A = 130 the SM basis is restricted 
to O57/2, 1^3/2, 1^5/2j2si/2 and Oh] 1 / 2 levels. The problem with these small basis sets is that the spin-orbit partners 
0/7/ 2 and 0<?7/ 2 have to be omitted [Tl|. The SM results then automatically violate the Ikeda Sum Rule (ISR) fljjj . 
while the QRPA satisfies it exactly. The Ikeda sum rule is: 

S- -S+ =3(N-Z), (10) 

A A 

5- = E<°< + iE(-) /i °--A'( fc )*+( fc )]E^(o*-(0]io i + >, (11) 

H k I 

For S+ the subscripts at the isospin rising and lowering operators are exchanged. 

The SM is in principle a better approach, if the single-particle basis is large enough. The ground state correlation 
in the QRPA are limited to the ring diagrams (see Fig.[I|, and thus to 'seniority' (the number of broken quasiparticle) 
0,4,8,12,... configurations, while the shell model can include seniority 0,4,6,8,10,12,... states. This has been 
discussed in publications of the Strasbourg-Madrid group [ll| . But the results [3, [ll[ show also that the philosophy 
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FIG. 2: (Color online) Contributions of the transforming neutron pair with different angular momenta J 7 ' to the total M°" 
calculated within the QRPA and different basis sizes for the 0vj3f3 decay 82 Se— 5- 82 Kr. The left bar is calculated with the same 
basis of four levels, IP3/2, 0/5/2, lpi/2 an d 0g$/2, used in the shell model calculations. The Ikeda Sum Rule (ISR) is exhausted 
by 50%. The second bar from the left includes in addition one, the I/7/2 level, of the two missing spin-orbit partners given 
for the 82 Se nucleus in ref. [t| for the shell model. The ISR is exhausted by 66%. The third bar from the left includes both 
missing spin-orbit partners O/7/2 and 0^7/2 amounting in total to 6 single-particle levels. The ISR is fulfilled by 100%. This 
leads to the increase in the neutrinoless matrix element from 1.12 to 4.07. The right bar represents the QRPA result with 9 
single-particle levels (I/7/2, 2p3/2, 1/5/2, %Pi/2, l<?9/2, 2^5/2, 3si/2, 2d 3 / 2 , lfl7/2-)- The matrix element gets only slightly increased 
to 4.27. The spin-orbit partners are essential to fulfill the Ikeda Sum Rule (ISR). In all four QRPA calculations the QRPA 
"renormalisation" factor g pp (given in the figure) of the particle-particle strength of the Bonn CD nucleon-nucleon interaction 
is adjusted to reproduce the experimental 2u/3/3 decay rates. 



of the RPA is approximately correct, since the ring diagrams give the most important ground state correlations. The 
contributions of seniority 6 and 10 to M 0v are suppressed compared to the others (see figure [4j. 

Figures [2] and [3] show the QRPA contributions of different angular momenta of the neutron pairs, which are changed 
in proton pairs with the same angular momenta. In figure [2] the left bar is the result for 82 Se obtained with the 
single-particle basis IP3/2, Of 5/2, lpi/2 an d Ogg/2 used in the SM. The ISR is exhausted by 50%. The second bar from 
the left represents the result with addition of the I/7/2 level. The ISR is exhausted by 66%. The third bar from the 
left shows the result obtained by inclusion of both spin-orbit partners O/7/2 and 0gg/ 2 missing in the four level basis. 
The ISR is 100% fulfilled. For the right bar the basis is increased to 9 single-particle levels for neutrons and protons 

(0/7/2, 1^3/2, O/5/2, lPl/2, 0# 9 /2, W5/2, 2s 1/2 , ld 3 / 2 , 0# 7 / 2 ). 

For 128 Te in figure [3] the left bar is calculated with the same five single-particle levels, 0<77/ 2 , 1^5/2, 2s 1 / 2 , 1^3/2 
and 0/in/2, used in the shell model calculations. The middle bar represents the results by inclusion in addition the 
missing spin-orbit partners Oc/9/2 and O/19/2, with 7 levels in total. The ISR is 100% fulfilled. This strongly increases 
the Qvfifl matrix element from 1.37 to 3.41. The right bar is the QRPA result with 13 single-particle levels including 
all the states from the N = 3 (p,f ) and the N = 4 (s,d,g) shells and the four additional levels 0hn/ 2 , 0/i 9 / 2 , 1/7/2 
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FIG. 3: (Color online) Contributions of different angular momenta associated with the transforming neutron pair for the QRPA 
with different basis sizes and with different angular momenta J 71 for the Ou/3/3 decay 128 Te— ^ 128 Xe. The left bar is calculated 
with the same five single-particle levels, 0177/2, W5/2, 2si/2, lc?3/2 and O/111/2, used in the shell model calculations. The middle 
bar includes in addition the missing spin-orbit partners Ogg/2 and 0/ig/ 2 , with 7 levels in total. This increases the neutrinoless 
matrix element from 1.37 to 3.41. The right bar is the QRPA result obtained with 13 single-particle levels including all states 
from the N = 3(p, /) and the N = 4(s, d, g) shells and the four additional levels O/in/2, O/1.9/2, 1/7/2 an d I/5/2 from the N = 5 
shell. The matrix element is only increased slightly to 3.82. The spin-orbit partners are essential to fulfill the Ikeda sum rule. 
The further increase of the basis from 7 to 9 levels produces only a small change from 3.41 to 3.82 for the 0uf5/3 transition 
matrix element. In all three calculations the QRPA "renormalisation" factor g pp (given in the figure) of the particle-particle 
strength of the Bonn CD nucleon-nucleon interaction is adjusted to reproduce the experimental 2u/3/3 decay rates. 



and I/5/2 from the N = 5 shell. The matrix element is increased only slightly to 3.82. The spin-orbit partners are 
essential to fulfill the Ikeda Sum Rule (ISR). The further increase of the basis from 6 to 9 levels produces a small 
change only from 3.41 to 3.82 for the 0^/3/3 transition matrix element. For all results in figures 2 and 3 the 2^/3/3 the 
QRPA "renormalisation" factor g pp of the particle-particle strength of the Bonn CD nucleon-nucleon interaction is 
adjusted to reproduce the experimental 2v[3/3 decay rates. The radius parameter is chosen to be ro = 1.2 fm. 

Figure 2] shows for 82 Se and for 128 Te the influence of different seniority states in the SM 8|. A seniority zero 
pair corresponds to a zero angular momentum neutron pair, which can be changed into a zero seniority proton pair. 
It yields a large positive contribution to the 0^/3/3 matrix element. The higher seniority states reduce the matrix 
elements as seen also in figures 2 and 3 for the QRPA. One remarks, that seniority 6 and seniority 10 states, which are 
not included in the QRPA, give only a very small or even zero contribution to the 0v[3/3 matrix elements in agreement 
with the RPA philosophy. 

The QRPA is restricted to the ring diagrams for the ground state correlations, while the SM takes into account 
all possible-many body states in a given single-particle basis. The drawback of the SM is the limited basis size. For 
the light double beta decay systems like 76 Ge and 82 Se the single-particle basis is restricted to the four levels [l-fUl 
IP3/27 O/5/2, lpi/2, 0(7g/2. (In ref. [9[ the I/7/2 level is also included for 82 Se.) Since the spin-orbit partners 0/ 7 / 2 , 0<7 7 / 2 
are missing, the Ikeda Sum Rule (ISR) is violated and one obtains only 50% (66% for five levels) of the the ISR (see 
table U). Including the spin-orbit partners into the basis O/7/2, IP3/2, 0/5/2, lpi/2, 0<7g/2, 0<7 7 /2 the ISR is fulfilled. It 
has been argued, that including the spin-orbit partner 0/ 7 / 2 ,0<7 7 /2 and moving the single-particle energy levels of 
these states to infinity also fulfills the ISR. This is correct, if one sums the states in the ISR to infinity. We have done 
a calculation with 6 levels moving O/7 /2 , 0<?7 /2 to 500 MeV. The ISR is fulfilled as expected, if one sums to energies 
larger than 500 MeV. However, the contributions to the ISR from low-lying states below about 10 MeV are within 1 
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FIG. 4: (Color online) The Ov/3/3 matrix element as a function of maximum seniority included in the SM wave functions for 
82 Se- 82 Kr and of 128 Te- 128 Xe. The contributions of seniority 6 and 10 states are appreciably smaller in agreement with the 
philosophy of the Random Phase Approximation. The data in the figure are taken from Caurier et al. ||. 

% the same as for the four level basis, and if one sums now the contributions to the ISR only from the physical states 
up to about 10 MeV, one obtains again only about 50 percent of the ISR. The nuclear structure properties of a six 
level basis with shifts of the two spin-orbit partners are markedly different from the calculation with the six levels 
with the correct single-particle energies. 

Table I shows the contributions from 'seniority' (the number of broken quasiparticles) zero and from 'seniority' 
different from zero to the total matrix element M 0l/ for different sizes of the single-particle basis in the QRPA 
calculations for 82 Se, 128 Te, 130 Te and for the large basis in 76 Ge, and for the small basis with 4 or 5 levels in the 
SM calculations [a-flll] for the two nuclei 82 Se, 128 Te. In all the QRPA calculations the "renormalisation factor" g pp 
of the particle-particle (pp) force multiplying the pp-nucleon-nucleon matrix elements of the G matrix of the Bonn 
CD interaction is adjusted to reproduce the experimental 2z//3/3 decay rates. In Ref. [l2[ the factor g pp has been 
adjusted together with the axial-vector charge qa simultaneously to the 2^/3/3 decay, to the single (3~ decay in the 
second leg and to electron capture in the first leg for three double beta decay systems. The results for the O^/3/3-decay 
matrix elements change with this overconstrained fit less than the general theoretical uncertainties (see figure [SJ. 
Only the Gamow- Teller transitions via the 1 + intermediate states depend sensitively on g pp , since the 1 + states in 
the intermediate nucleus are close to a "phase transition" to "static Gamow- Teller deformations". As a function 
of the other particle-particle multipole parts of the force the results change only very weakly. Therefore, possible 
uncertainties of g pp , which affect a small contribution of the 1 + multipole, do not affect appreciably the total 0^/3/3 
matrix element. 

Table II shows the percentages of the different 'seniority' (the number of broken quasiparticles) admixtures in the 
QRPA ground states obtained with the Bonn CD NN force for different nuclei and the basis with 9 and 13 single- 
particle levels defined in the caption of table I. The formalism for the determination of the relative percentages for the 
different seniorities is described in ref. [l3j]. Making use of eqs. (J3j) , (J4J) and (0, the admixtures of different 'seniorities' 
(the number of broken quasiparticles) can be determined by the following eqs. (|12[) to (|14p . The excitation of + 
unbroken quasiparticle pairs should be excluded in the sums of eq. (fl~3"|) . 



7 



TABLE I: Contribution from seniority zero and seniorities different from zero to the total value of the Ov(3(3 transi- 
tion matrix elements obtained for different nuclei in the QRPA and the SM with different basis sizes. The last column 
gives the exhaustion of the Ikeda Sum Rule (ISR) in percents (see eqs. (|10[1 . ([TT]) ) . The basis sets used are: 4 levels: 
1^3/2, 0/5/2, lPi/2,0ff 9 /2 ; 5 levels: 0g 7 / 2 , W5/2, 2s 1/2 , ld 3/2 , 0h 11/2 ; 6 levels: 0/ 7/2 , lp 3 / 2 , 0/ 5 / 2 , lpi/2, 0g g/2 , 0g 7/2 ; 7 levels: 
O37/2, ld B/2 ,2si/ 2 , ld 3/2 ,0h 11/2 ,0g 9/2 ,0h 9/2 ; 9 levels: 0/ 7/2 , lp 3 /2,0/ 5 /2, lfi/2, 0cjg/ 2 , ld 5/2 ,2s 1/2 , ld 3/2 ,0g 7/2 ;13 levels: All lev- 
els of N = 3 (p, /) and the N = 4 (s, d, g) shells and 0/in/ 2 , 0/ig/ 2 , 1/7/2 and 1/5/2- Four levels (for: 76 Ge, 82 Se) and five levels 
(for 100 Mo, 128 Te, 130 Te) correspond to the SM basis sets [8Blll]. Since the level four and five basis sets do not contain all spin 
orbit partners, the ISR is strongly violated and only exhausted by 50 and 60 percents for these sets, respectively. 



Nucleus 


Basis Size 


Method 


seniority: s = 


seniority: s / 


total 


Ikeda SR percents 


76 Ge 


9 


QRPA 


12.5 


-7.5 


5.0 


100 


82 Se 


4 


QRPA 


6.6 


-5.5 


1.1 


50 


82 Se 


4 


SM 


7.8 


-5.8 


2.0 




82 Se 


5 


QRPA 


8.9 


-6.0 


2.9 


66 


82 Se 


5 


SM 


+ /7/2 




2.5 




82 Se 


6 


QRPA 


10.7 


-6.6 


4.1 


100 


82 Se 


9 


QRPA 


11.8 


-7.5 


4.3 


100 


100 Mo 


13 


QRPA 


16.3 


-12.6 


3.7 


100 


i28 Te 


5 


QRPA 


9.7 


-8.3 


1.4 


60 


128 Te 


5 


SM 


10.6 


-8.4 


2.2 




128 Te 


7 


QRPA 


13.7 


-10.3 


3.4 


100 


128 Te 


13 


QRPA 


16.8 


-13.0 


3.8 


100 


130 Te 


5 


QRPA 


8.8 


-7.5 


1.3 


60 


130 Te 


7 


QRPA 


12.2 


-9.0 


3.2 


100 


130 Te 


13 


QRPA 


14.9 


-11.1 


3.8 


100 



TABLE II: Contributions (in %) of different 'seniorities' (the number of broken quasiparticles) to the normalization of the 
QRPA ground states for different nuclei. The 9 and 13 level single-particle basis sets are defined in the caption of table I. The 
renormalisation factor g pp of the Bonn CD NN force is adjusted to reproduce the experimental 2^/3/3 half lives. A larger basis 
and also a stronger NN force produces larger admixtures of higher 'seniorities'. 



nucleus 


basis 


s = 


s = 4 


s = 8 


s = 12 


s = 16 


V0 Ge 


9 


55.0 


32.9 


9.8 


2.0 


0.3 


76 Se 


9 


58.7 


31.3 


8.3 


1.5 


0.2 


82 Se 


9 


56.4 


32.3 


9.3 


1.8 


0.3 


82 Kr 


9 


53.5 


33.5 


10.5 


2.2 


0.3 


128 Te 


13 


51.5 


34.2 


11.3 


2.5 


0.4 


128 Xe 


13 


39.9 


36.7 


16.9 


5.2 


1.2 



^T = E^ X ™> o^ = 5>™(A7)-\ (12) 

m 

w ( s = 4 ) = iEi c '^i 2 ( 13 ) 

a/3 

w(s = 8) = ^w(s = 4) 2 ; w(s = 12) = ^w{s = 4) 3 ; w(s = 16) = ^w(s = 4) 4 (14) 

A comparison of the QRPA and the SM results for the Gamow Teller contributions of the different angular 
momentum pairs to the transition matrix elements is shown in figure [5] The basis has been chosen for the QRPA 
to be the same small basis as for the Shell Model. ( For S2 Se : 1^3/2, OI/5/2, 1.P1/2, O99/2 an d for 130 Te : for protons 
O59/2, 057/2, 1^5/2, 14/2, 2si/ 2 and for neutrons 0g 7 / 2 , W5/2, W3/2, 2si/ 2 , 0h n / 2 -) 

In the last ten years P. K. Rath and coworkers [TB4I3 nave published a whole series of papers (see references in 
Ref. [3|) on 2z//3/3 decay and, since 2008, also on 0v/3f3 decay, in which they used a simple pairing plus quadrupole many 
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body Hamiltonian of the Kumar and Baranger type 24[ to calculate the neutrinoless double beta decay transition 
matrix elements with angular momentum projection from a Hartree-Fock-Bogoliubov (HFB) wave function after 
variation. Schmid [2l[ did show, that with the assumption of a real Bogoliubov transformation (real coefficient A 
and B) 



M 

a£ = ^ A ta c\ + B ia a (15) 

i=l 

and no parity mixing, only + , 2 + , 4 + , . . . . nucleon pairs and excited states are allowed. In this approach one can 
describe rotational bands and can change neutron pairs with angular momentum + , 2 + ,4 + , . . . into a correspond- 
ing proton pair with the same angular momentum and parity in neutrinoless double beta decay for the ground to 
ground state transition in even mass nuclei. The proof uses the Bloch-Messiah [25| decomposition of the Bogoliubov 
transformation. After the first Bloch-Messiah transformation 

M 

K = Y. D ^l (16) 

i 

one can write the HFB state as a BCS wave function: 

\HFB)= \{ [u a + v a bibl}\0) (17) 

a,m a >0 

The bar indicates the time reversed nucleon state. The neutron pairs to be changed into two protons in a + — > + 
Qvj3f3 transition can be written for axial-symmetric deformed nuclei: 

T=p,n i<k J 

■MDl a D k , a )[l + (-)<«+<*+'] + i$(D* a D k , a )[l (-) k+h+J ]][clcl}{ : l (18) 
J2 E(V2, 1/2, T\ - 1/2, 1/2, )(-)^™ (Ji, 3k, J\ma, -m a , 0) 

i,k J,T 

■[M(D* a D ka )[l + (-y^ k+J] + i^ {DLDka)[1 _ (-)'«+ l *+ J ]][cJ c t]g. 

If only real coefficients are admitted and no parity mixing is allowed in the Bogoliubov transformation (| 15[) (l4l-fl6| , 
the expressions for the neutron pairs to be transformed into a proton pair for a + — > + transition reduces to: 

blbt = £[1 + S i!k ]- 1 J2(-y k+jk - m "(ji,jk, J\m a ,-ma,0) ■ [(D* a D ka )[l + {-) J ]][c\cl]{ fl 2T . (19) 

i<k J 

The neutron pairs, which can be transformed into a proton pair for + — > + transitions are restricted to angular 
momenta + , 2 + , 4 + , 6 + , . . . . 

The QRPA and the SM do not have this restriction. Figure [6] shows the contributions of the different angular 
momentum and parity pairs to the 0i//3/3 matrix elements in 76 Ge, 100 Mo and 130 Te calculated within the QRPA. 

Figure [7] shows on the other side, that the projected HFB approach is restricted to contributions of neutron pairs 
with angular momenta + , 2 + , 4 + , . . . . In addition, one sees that the contributions of transition of higher angular 
momentum neutron to proton pairs 2 + ,4 + ,... are drastically reduced compared to the QRPA and the SM see 
fig. [7] The reason for this is obvious: in a spherical nucleus the HFB solution contains only seniority zero and no 
stronger higher angular momentum pairs. The double beta decay system considered has only a small deformation 
and thus a projected HFB state is not able to describe an appreciable admixture of higher angular momentum pairs 
for + — > + transitions as can be seen from the wave function (fl7|) after the first Bloch-Messiah transformation 
(25| . The higher angular momentum contributions increase drastically with increasing intrinsic quadrupole and 
hexadecapole deformations of the HFB solution. 

The results in figure [Jjare calculated by K.W. Schmid [26| within the HFB with the angular momentum projection 



before variation with an improved Gogny force 27J adjusted in a global fit to properties of many nuclei. 

To have also l + ,3 + ,5 + , . . . neutron pairs contributing one has to use a Bogoliubov transformation with complex 
coefficients A and B [see equations (fT5")) and (|19l) ]. To have also CP, 1~, 2 _ , 3— , 4 _ , 5 _ , . . . one has to allow parity 
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mixing in the Bogoliubov transformation ( I15[) . But even allowing all different types of angular momentum and parity 
pairs one would still have an unnatural suppression of the higher angular momenta especially for smaller deformations. 
This handicap could probably be overcome by a multi-configuration HFB wave function [2l[ with complex coefficients 
and parity mixing in the Bogoliubov (|15l) transformation. 

The IBM (Interacting Boson Model) [18| can only change + (S) and 2 + (D) fermionic pairs from two neutrons into 
two protons. In the bosonization to higher orders this leads to the creation and annihilation of up to three "s" and "d" 
boson annihilation and creation operators in Ref. [l8j . But all these terms of equation (18) of reference 18] originate 
from the annihilation of a + (S) or a 2 + (D) neutron pair into a corresponding proton pair with the same angular 
momentum. The higher boson terms try only to fulfill the Fermi commutation relations of the original nucleon pairs 
up to third order. The IBM can therefore change only a + or a 2 + neutron pair into a corresponding proton pair. 

III. CONCLUSIONS 

The Shell Model (SM) @, Hoi . ITT1 ] is in principle the best method to calculate the nuclear matrix elements for neutri- 
noless double beta decay. But due to the restricted single-particle basis it has a severe handicap. The matrix elements 
in the 76 Ge region are by a factor 2 smaller than the results of the Quasiparticle Random Phase Approximation 
(QRPA) [Ml. 

With the same restricted basis as used by the SM the QRPA obtains roughly the same results as 
the SM (see figure 5), but the Ikeda sum rule [T^ j gets strongly violated due to the spin-orbit partners missing in 
the SM single-particle basis. Future increases in computer performance and perhaps also more powerful methods of 
programming might bring improvements for the SM calculations of M 0u . 

The SM can admix seniority 0, 4, 6, 8, 10, . . . configurations to the ground state, while the QRPA is restricted in the 
ground state to the "seniority" (the number of broken quasiparticles) 0, 4, 8, 12, . . . configurations (see figures [T]and 
But the SM results show also, that the seniority 6 and 10 configurations, which arc not included in the QRPA 
ground state, contribute just a little to the total neutrinoless double beta deca y m atrix element (see figure [4]). 

The angular momentum projected Hartee-Fock-Bogoliubov (HFB) method [15[ is restricted in its scope. With a 
real Bogoliubov transformation without parity mixing (1151) one can only describe neutron pairs with angular momenta 
and parity + , 2 + , 4 + , 6 + , . . . changing into two protons for ground state-to-ground state transitions. The restriction 
for the Interacting Boson Model (IBM) [l8| is even more severe: one is restricted to + and 2 + neutron pairs changing 
into two protons. 

A comparison of the 0v/3[3 transition matrix elements calculated recently in the different many body methods: 
QRPA, SM, HFB and IBM is shown in Fig. M 
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FIG. 5: Gamow Teller (involving the spin matrix it) contributions of different angular momenta associated with the transforming 
neutron-pair for the QRPA and the SM in 82 Se (upper panel) and 130 Te (lower panel). Since the Gamow Teller part is shown 
here only, the results are slightly different from figure [2] and table [I] where the total matrix element is given. The SM results 
shown in this figure do not include the effects of the higher order currents [1, QlJ, while the QRPA ones do include them. The 
basis is the same small Shell Model basis for QRPA and the SM (four levels for 82 Se and five for 130 Te). 
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FIG. 6: (Color online) Contribution of different angular momenta for both parities (positive above and negative below) of the 
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FIG. 7: Contributions of neutron pairs with different angular momenta to the neutrinoless double beta decay transition matrix 
elements for 76 Ge — >- 76 Se calculated from a HFB wave function with angular momentum and proton and neutron particle 
number projection before variation. The Fermi, the Gamow- Teller and the total contribution including the tensor part as 
defined in eq. ([8]) are separately given. The nucleon-nucleon interaction is an improved Gogny type force [27|. The deformations 
Pae = —0.08 and j3s e = 0.11 correspond to the minima of the projected HFB total energy. The results are qualitatively and 
almost quantitatively the same for the experimental deformation from the Coulomb reorientation effect: /3ge = 0.16, /3s e = 0.10 
and also for different forces. The angular momenta of the neutron pairs are restricted to + , 2 + , 4 + , . . . . In addition the 
contributions of higher angular momentum neutron pairs 2 + , 4 + , . . . are drastically reduced compared to the QRPA and the 
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FIG. 8: (Color online) Neutrinoless double beta decay transition matrix elements for the different approaches: QRPA @, @|, 
the SM [EES, El, the projected HFB method QjJ and the IBM [lq ]. The error bars for the QRPA are calculated as the highest 
and the lowest values for three different single-particle basis sets, two different axial charges gA = 1-25 and the quenched value 
qa = 1-00 and two different treatments of short range correlations (Jastrow-like [28[ and the Unitary Correlator Operator 
Method (UCOM) [2S|). The radius parameter is as in this whole work ro = 1.2 fm. 



